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Abstract 

We examine the one-dimensional motion of two similarly charged particles under 
the influence of only two forces, i.e. their Coulombic repulsion and their gravitational 
attraction, using the relativistic equation of motion. We find that when the rest mass 
of the two particles is sufficiently small (~ 0.4 eV/c 2 ) and the initial Coulombic po- 
tential energy is sufficiently high (~ m p c 2 , where m p is the proton mass), then the 
strong gravitational attraction resulting from the relativistic particle velocities suffices 
to counterbalance the Coulombic repulsion and to cause stable periodic motion of the 
two particles. The creation of this confined oscillatory state, with a rest mass equal to 
that of a proton, is shown to be consistent with quantum mechanics by examining the 
particle dc Broglic wavelength and the Klein-Gordon and Schrodinger equations. It is 
shown that the gravitational constant can be expressed in terms of the proton mass and 
charge, the vacuum dielectric constant, the Planck constant and the speed of light. It is 
also shown that gravity can cause confinement of light (~ 0.1 eV/c 2 ) neutral particles 
(neutrinos), or pairs of a neutral and a charged light particle, in circular orbits of size 
0.9 fm and period 10~ 24 s forming bound neutral or charged hadron states. 
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1 Introduction 



The possibility that gravity may have a significant role at short , femtometer or subfem- 
tometer distances has attracted significant interest for years and the 

potential role of special relativity [8| as well as the feasibility of developing a unified gauge 
theory of gravitational and strong forces 01 have been discussed. 

Gravitational forces between small particles at rest are entirely negligible in comparison 
with Coulombic forces. Thus the gravitational attraction between two protons at rest is 36 
orders of magnitude smaller than their Coulombic repulsion, i.e. Gm 2 = 8.03 • 10~ 37 (e 2 /e) 
where we denote e = 47re D . 

Due to this enormous 36 orders of magnitude gap, little attention has been focused 
on gravitational forces between fast moving particles with velocities very close to c. For 
the laboratory observer in frame S (Figure 1) the one-dimensional relativistic equation of 



motion of a particle with rest mass m is [9|; [lCj] : 



dp d{^m v) 3 dv dv 
dt dt ^ ° dt ^ dt 

where mp = 7 3 m Q is the longitudinal mass of the particle, 7(= (1 — v 2 /c 2 ) -1 / 2 ) is the Lorentz 
factor and v =v is the velocity of the particle relative to the laboratory observer 



BEloi 



As 

([T]) shows, mpjs the ratio of force divided by acceleration, thus it is the inertial mass of the 
particle |g;|9|;ll0l|. i.e. the quantity defined and measured as the mass of all bodies subject to 
gravity [9cllCl|. It is on the basis of this quantity that Newton's r~ 2 law has been formulated 
and Newton's constant has been measured. In view of the well proven equivalence principle 
it is thus equal to the gravitational mass in Newton's r~ 2 gravitational law [8]. It is also 
the only mass value for the moving particle B which the laboratory observer S can judge on 
the basis of his force and acceleration observations, the rest mass m of the moving particle 
cannot be measured in frame S, only m.£ is measurable in the laboratory frame. Thus it is 
this quantity, rather than the rest mass m , which is available to the laboratory observer 
to use in Newton's gravitational law applied to the one-dimensional motion of the moving 
particle, i.e. 



F G = ^ = m^ = ^ = ^f- ( 2 ) 



dp dv Gm 2 Gm 2 '-/ 6 
dt dt r 2 r 2 

For particles with non-relativistic velocities this subtle difference is negligible. Also at 
a first glance the 7 6 term appears insufficient to make Fq significant relative to Coulombic 
forces unless the velocity v is high enough to bring the Lorentz factor 7 close to 10 6 . However 
such relativistic velocities can in principle be easily reached by small particles under the 
influence of their Coulombic attraction or repulsion. 

Thus in this work we examine the seemingly very simple problem of the one-dimensional 
motion of two similarly charged small particles under the influence of only two forces, i.e. 
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Figure 1: Schematic of the two charged particles with rest mass m , of the symmetry axis 
and of the two reference frames S and S'. Inset summarizes the key 0; OjJ relativistic 
relationships between the observations in frames S and S 1 ) 7 3 (r)m is the longitudinal mass 
[lo| at displacement r; The same force value, F, is observed in both frames, but the two 
observed accelerations, a' and a, have a ratio of 7 3 0; 10]. 



their Coulombic repulsion and gravitational attraction using the relativistic equation of 
motion and the Coulomb and Newton laws in order to examine under what conditions the 
two particles will not simply escape from each other, but may collide with each other or 
form a stable oscillatory state. At the end we show the consistency of the results with 
quantum mechanics by examining the de Broglie wavelengths of the two particles and the 
Klein-Gordon equation of quantum mechanics which allows for the use of the Schrodinger 
equation under relativistic conditions. 

We focus on examining the conditions for which the rest mass of the oscillatory two- 
particle system corresponds to that of a hadron such as a proton. There has been some 
preliminary work in this area treating hadrons as standing waves 11] or strings [l^] and 
leading to the same analytical expression for G derived here, but the exact origin of the 
creation of the standing wave bound state was unclear. The formation of hadrons via 
condensation of smaller 
in the QCD theory 131 : 



particles, i.e. of the gluon-quark plasma, is commonly analyzed 

3; 0; 0; 0; 

This condensation occurs at the transition 



temperature of QCD which is given by Tq = 151(6) MeV in the kT scale 17J, i.e. it 
corresponds to a particle energy of approximately 150 MeV 17]. 



2 Reference frames 

As is common practice 

HQ 

, we examine the motion of each particle (e.g. particle B) using 
two reference frames (Figure 1). The laboratory frame S is at zero velocity with respect to 
the center of mass of the two particles, while the instantaneous rest frame S' has a velocity 
with respect from S equal to that of particle B jfll: [lo|. 

The inset of Figure 1 summarizes the basic elements of special relativity [s; 10] regarding 
the values of force, velocity and acceleration of particle B observed in the two reference 
frames. The only one of these basic elements needed for the present analysis is that, although 
the force value, F, observed in both reference frames is the same 0; [h]], the acceleration, 
a, measured by the laboratory observer S is 7 3 times smaller than the acceleration, a', 
measured in frame S' and thus the mass value judged by the laboratory observer is 
7 3 times larger than the mass m judged in frame S', which is the rest mass m of the 
particle since the particle is at rest (u' = 0) in frame S' 0; 0]. Thus the longitudinal 
mass mi = 7 m is the only mass value observed by the laboratory observer and thus, in 
view of the equivalence principle 0], the only one available to use in Newton's gravitational 
law for the one-dimensional particle motion, as already discussed. The rest mass m of the 
accelerating particle is not measurable by the laboratory observer. This simple but key 
point is worth emphasizing because what follows after equation ([5]) in treating the stated 
one-dimensional particle motion problem is then basically simple energy and momentum 
conservation and the corresponding simple algebra. 



3 Initial conditions 

We denote r Q the initial particle distance and we define x = r/r a . Since we are interested in 
the possible formation of a hadron with charge e, i.e. a proton, we consider due to charge 
conservation, that the charges of the two particles of rest mass m each are q\ and q^ with 
Qi+Q.2 = e- Denoting qi/e = Q\ and 02/e = Q2 and Q1Q2 = Q (< 1/4) it follows that the 
initial Coulombic potential, Vc, , is given by: 

Vc,o = ^ (3) 
er Q 

We denote: 

Vn 



m p c z 

where m p is the proton mass, thus: 



(4) 



to = ~^- 2 = (£) \ q = (?) \ c (a/2«) (5) 
pem p c 2 \2p J \p J 

where A c (= h/m p c = 1.32 fm) is the proton Compton length, a(= e 2 /ech = 1/137.035) 
is the fine structure constant (so that a/2-7r = e 2 /ech), and the length X q (= 2e 2 / sm p c 2 = 
2A c (a/27r) = 3.06 • 10~ 3 fm) is of the size range of quarks (thus we use the subscript q) and 
is twice what is commonly termed "classical radius" of the proton. 
Since x = r/r Q it follows from © that: 

-(?)(*) 

Thus the parameter p{= Vc, /mpC 2 ) defines via ([5]) the initial (t = 0) particle distance 
r Q and thus the initial dimensionless distance x Q (= 1). 

In choosing the initial velocity, v Q , or equivalently the initial kinetic energy K Q {= — 
l)m c 2 ) of each particle, we have examined two general cases: (a) 2K Q = (1 — p)m p c 2 , so 
that 2K Q + V c o — TfipC , which we denote case A. (b) K Q — 0, thus 7 Q — 1, which we denote 
case B. 

The two cases coincide for p = 1 . 

4 Energy conservation 

If due to the action of gravity the two particles form a bound state then the total energy 
2E(x) = 2j(x)m c 2 of the two initial particles becomes the rest energy, R(x), of the confined 
state. In this case the two particles can be viewed as partons of the confined hadron state 
in Feynman's parton model or, as shown later, as quarks and gluons in K. Johnson's bag 
model or in the standard model. If the rest mass R(x) is at some point x p equal to the rest 
energy of a proton, i.e. R(x p ) = m p c 2 , then denoting j(x p ) = j p it is: 

m p = 27 p ra ; >y p = m p /2m (7) 

It is worth noting that since in general the total particle energy E{x) is related to its 
momentum p(x) via: 

E 2 {x) =m 2 c 4 +p 2 (x)c 2 (8) 
and since R(x) = 2E(x) it follows: 

R 2 (x) = Am 2 c 4 + Ap 2 (x)c 2 (9) 



We start the derivation by using energy conservation: 



R(x) = 2-y(x)m c 2 = 2m c 2 + 2 (7(2;) - 1) m c 2 = 2m c 2 + 2K(x) 
= 2m c 2 + 2K + V c ,o + V G ,o ~ V G (x) - V G (x) 



(10) 



where K Q is the initial kinetic energy of each particle, V G and V G)Q are the Coulombic and 
gravitational potential energies at x = 1 (i.e. at r = r a ), V G (x)(= V G , /x) and V G {x) are 
the Coulombic and gravitational potential energies at x and i^(x)is the kinetic energy of 
each particle at x. 

Choosing to analyze first the case of initial conditions A, i.e. 2K Q = (1 — p)m p c 2 so 
that Vc : o + 2K a = pm p c 2 + (1 — p)m p c 2 = m p c 2 and dividing (jTUJ) by 2 r y p m c 2 = m p c 2 one 
obtains: 

i{ x )h P = K x ) = 1 - p (~ -y( x )J (11) 

where we have accounted for V G , + K a = m p c 2 , for Vc(x) = V GjQ /x and for m <C m p , 
i.e. for 7 P S> 1 which is the case as shown later, and where we have defined y(x) = 
-{V G (x) - Vg,o)/Vc, and p{x) = j(x)/"/ p = R(r)/m p c 2 . 

Note that for zero relative velocity of the two particles it is y{x) = £/x, where £ = 
eGm 2 /e 2 « 10~ 36 , but for relativistic velocities the function y(x) becomes significant and 
has to be determined (section 4.1). 

It follows from (jlip that when p{x) = 1 or, equivalently, R(x) = m p c 2 then: 

y(x p ) = l/x p ; -(y G (x)-V Gt0 ) = V c (x p ) (12) 
In case B equation (fl~Qj) becomes: 



R(x) = 2>y(x)m c 2 = 2m a c 2 + 2(7(2;) - l)m c 2 = 2m c 2 + 2K(x) 

= 2m c 2 + V c ,o + V GtQ - V c (x) - V G (x) (13) 

and accounting for — V Gj0 « 10~ 36 and dividing by 2^ p m c 2 = m p c 2 one obtains: 



l( x )/l P = P( x ) = P fl - ^ + Vg(x)/Vc, 
7" 1 + p (l - ~ + y(x)J =p(l-~+ y{x)J 



(14) 



where the last equality holds since, as already noted, m <C m p , thus 7 P S> l.For p = 1 
equations (fTTj) and (fl~4"|) coincide. 



4.1 The variation of gravitational potential energy with distance 

The gravitational potential energy, V G (x), can be computed using ([2]) for any x from: 

Gm 2 rx ' =x 



- (VgO) - V G ,o) 



-F G {x') d(r x') 



x'=l 



x'=l 



J"(X_ 
r l2 



■dx' 



(15) 



where x' denotes the dummy variable. 
Using the definition of y(x), it follows: 



V{x) = -(V G (x)-V Gi0 )/Vc, = -(V G (x)-V Gy0 )/Q(e 2 /er ) = -(V G (x)-V Gy0 )/pm p c 2 (16) 



and thus: 



y(x) 



(V G (x) - V G ,o) Gm 2 



x'=x ^Q(„i\ 

7 [X) dx> 



Q(e 2 /er ) Q(e 2 /e) J x , =1 x' 2 
and using m p = 2m ^ p one obtains: 

eGml 2 /-'=- 7 6( x ') , sGm 2 p 4 7 6 (^)/7 P % , 



AQe 2 Ip J x , =1 x' 2 AQe 
Using (fTTj) to express 7(2/) in (fT8j) one obtains: 



x'=l 



x 



12 



(17) 



(18) 



with: 
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dx' = 46 



x'=l 



x'=l x ' 

[l-pjl/x'-yjx'))? 



■dx' 



46 - £Gm '~ 4 



4Qe 



% or 46 = £ 
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(19) 



(20) 



The third equality (|20p defines £, which is the constant we aspire to determine, while 
the last equality shows that the parameter b is uniquely defined by f; and by the choice of 
the mass and charge of the two particles. 



The measurement of G is commonly carried out using metal rod torsion balances 19|; 
20j; [2l( On the basis of the CODATA recommended [19J] experimental values of e, G, m p 
and e (Table 1) it is: 



£ exp = 8.09335 • 10~ 37 = 0.134357(a/27r) 



12 



(21) 



Thus if, as an example, one uses this experimental value for £ and also chooses Q = 
2/9 (corresponding to q\je = 2/3 and qije = 1/3, the charges of u and d quarks) and 
m = (l/4)(a/27r) 3 m p = 0.3675 eV/c 2 , (close to the estimated heaviest neutrino mass of 
0.4 eV/c 2 22|) then (|20p gives b = 0.6046, which lies in the range of b values leading to 
particle confinement as shown below. 

Upon differentiation of (|19|) using Leibnitz's rule one obtains: 

^ = 4fr (w/x t Mx))6 ; with B - c - = ; case A ( 22 ) 

ax x l 
In case B in view of (|14p the differential equation (1220 becomes: 



d^x)_ = 6 (J_J^x + yf wUhBC y(1) = . case B (23 ) 
ax x z 

When p = 1 equations (f2"2"j) and (|2"3"j) are identical. 

Once 6 and /O have been chosen and the differential equation (|22p or (|23p has been solved 
numerically to obtain y(x), then recalling in case A or (|14|) in case B and the definition 
of y(x) it is: 

fj,(x) = 1 - p f- - y{x)^j ; case A (24) 

fj,(x) = p ( 1 \-y(x)) ; case B (25) 

V x J 

and thus in either case the function p{x) is determined. Thus both the potential energy 
profiles and also the force profiles are readily determined. Thus from Coulomb's law the 
Coulombic potential energy is: 

V C (x) = V c ,o/x (26) 
and from the definition of y(x) = —(Vg(x)—Vg,o)/Vc,o t ne gravitational potential energy 

is: 

- (Vg(z) - V G ,o) = VcM*) (27) 

Examples of the p(x), Vcix) and Vg(x) profiles for p = 1 (thus K Q = and —Vg,o = 
£Vc,o ~ lQ~ S6 Vc,o = 0) an d two different b values (b = 2/9 and b = 8/15) are presented in 
Figures 2 and 3. These two b values have been chosen to lie below and above a critical b 
value (b c = 0.32 = 8/25) above which, as shown below, particle confinement occurs. 

As shown in Figure 2a for 6 = 2/9 the gravitational potential Vq energy approaches a 
limiting value of —220 MeV for x — > oo while for b = 8/15 (Fig. 2b) — Vn(x) approaches in- 



finity at a distance of 5 • 10 f m (asymptotic freedom behaviour LJ; Il5l;ll6|]). This bevavior 



is typical for b values above b c = 8/25 and leads, as shown below, to particle confinement. 




Figure 2: Coulombic (Vc) and gravitational (Vq) potential energy profiles for p = 1, from 
flUD and (J2ZD, with 6 = 2/9 (a) and 6 = 8/15 (b). 



Figures 3a and 3b show the corresponding p,{x) and R(x) profiles for b = 2/9 and b = 8/15 
respectively. One observes that in the former case R(x)/c 2 approaches an asymptotic value 
(1157 MeV/c 2 ), interestingly very close to the masses (1115 to 1197 MeV/c 2 ) of the £, E + , 
S° and S~ spin 1/2 baryons [22j]. In the latter case (b = 8/15), fx(x) and R(x) rise steeply 
after x p (where p,(x) = 1 and R(x) = 938.272 MeV) at some point xt which, as shown 
below, corresponds to the terminal (maximum) distance during the oscillations of the two 
particle system. It is again interesting to note (Fig. 3b) that the rest masses of all the spin 
3/2 baryons (Q~ , H*, S* and A) [22] all correspond to x values between x p and xt- 
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Figure 3: Rest mass R(x) and corresponding /i(x) profiles for p = 1, from (|24p or (|25p . with 
6 = 2/9 (a) and b = 8/15 (b). Comparison with the rest masses of the spin (1/2) (a) and 
the spin (3/2) baryons 



4.2 Force profiles, maximum momentum and confinement conditions 

The force profiles Fc(x) and Fq{x) can be obtained either via differentiation of the potential 
energy profiles (|26p and (|27p or directly from Coulomb's and Newton's Laws, i.e.: 



Qe 2 Qe 2 V c , 2p 2 m p c 

Fc(x) = j = 2 = 2 = 7T\ 2 V 28 ) 

where in the last equality we have used ([3]), ([!]) and ©, and thus also using ([20]) : 



Gmg 7 |(x) Gmg 7 6 (x ) Gm^ (Qe 2 / £ ) 7 6 (x) 

■'Gl 1 ) = ~3. = 3 " 



r 2 r 2 47 2 (Qe 2 /e) r r x 2 

= 4fe 7 -6 ^o7 6 (^) = 46 ^6 (x)Fc(x) (29) 

Therefore the net force i^x) is given by: 

F(x) = F c {x) + F G {x) = (1 - ^(x))^^) (30) 

and thus a necessary, but not sufficient, condition for particle confinement is that there 
exists a distance x m , such that: 

^(x m ) = 1/46 (31) 

At such a point it follows from (pQ) and (I30p that the momentum, j?(x), exhibits a local 
maximum. 

Although not directly needed here, we note that numerical integration of (|22p . i.e. case 
A, choosing various p and b values shows that for all p values above 0.7 the condition (|3ip 
is equivalent to: 

b > 6 min = 0.136244 (32) 
while in case B the maximum momentum condition (|32p is replaced by: 

bp 6 > 6 min = 0.136244 (33) 

The force profiles for p = 1 and 6 = 2/9 and 8/15 are given in Figure 4. One observes that 
in the latter case, asymptotic freedom behavior is obtained [l^: fla: [l^| . i.e. the attractive 
force —Fq, which is small at short distances, approaches infinity with increasing distance. 
Thus in this case (p = 1,6 = 8/15) the two particles cannot escape from each other. 
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Figure 4: Coulombic (Fc), gravitational (—Fg) and total (F = Fq + Fq) force profiles, 
from eqs. ([28]) to ([30]) . for b = 2/9 (a) and b = 8/15 (b). In the latter case, asymptotic 
freedom behaviour [lj; EEI; 03] is obtained. 



5 Momentum profiles and sufficient confinement condition 



Since the force expressions (]28p to (|30p are valid in both cases A and B, the equation of 
motion and momentum profiles discussed in this section are the same for both sets of initial 
conditions A and B. The only difference is that in case B the function p,(x) has to be 
computed from (j25H rather than from (|24p . 

A necessary and sufficient condition for particle confinement is obtained by using ([2]) 
and (I30p and examining the equation of motion of particle B: 



^ = F(x) = Fc(x) + F G (x) = F c (x) [l - Ab^(x)] 
where p(x) = ^(x)m v(x). This equation which can also be written as: 



(34) 



dp dr , dp(x) . , „. . dp(x) ±\ a Fn(x) , , n 

/ ' T+ = F x ; -^V"^) = F 0*0 ; -^r 1 = — ; x i/2 1 - * ( 35 ) 

dr dt r Q dx dx n -2V' 



where we have used ([5]) with Q = 1/4 (identical particles), thus, r Q = (l/8p)X q , and 
have also used the definition of the Lorentz factor , 7 i.e. (v/c) = ±(1 — 7-" 2 ) 1 / 2 . Using 
([28]) to express Fc(x) one obtains: 

Integration of this equation gives the momentum profile p(x), i.e.: 




Figure 5: Particle momentum profiles from (138(1 for p = 1 and various values of b and 
corresponding rest particle mass Tfi Q) also showing the definitions of x m and xt for the case 
b = 0.4. The momentum exhibits a local maximum and minimum for b > b m [ n (m < 
"T-o,max = 0.516 eV/c 2 ) and the particles are confined (finite xt) for b > b c = 8/25 (m < 
mo, c = 0.418 eV/c 2 ). It is b = (eGm 2 /16Qe 2 )(m p /2m ) 4 



, / m f 1-46/iV , , 

p{x)=p + pm p c - dx (37) 

Jx'=i (1 --f- 2 (x')) 11 V 2 

where p D is the momentum at x = 1. This equals zero in case B (X G = 0) and can be 
easily shown from ([9]) to equal, to a very good approximation, m p c(l — p) in case A. Thus 
using also (fTHl) to express 7 P (= (46/^) 1 / 4 ) one obtains: 



f x=x l-46u 6 (x') 
2(1 - p) + 2p / p \ dx' ; case A 

Jx'=l (l- 7 - 2 (x0) 7 X' 2 



{m p /2)c ' J x , =1 (l- 7 -2( x /)) 

|p(x)| /-^-^ 1 -46/i 6 (x') , , „ 
! ' ! - 2p / — \fz dx' ; case 5 (38) 



(m p /2)c J x , =l (l- 7 -2( x /))i/2 x /2 

Thus, focusing on case B, a necessary and sufficient condition for particle confinement 
and establishment of a self-sustained oscillation is that there exists some finite positive 
distance xt at which p(xt) = 0, i.e. the momentum vanishes. At this point it is dp(x)/dt < 
and thus at xt the velocity becomes negative until x = 1 is reached thus completing the 
cycle (Figure 5). 

In this case the period of the oscillation, T OSCl can be computed using ([54^1 . i.e. 



e.g. for case B, using ([35]) : 

^ rx'=x T ^ x > 
T0Sc/2= 8icJ xt=1 (l- 7 -2 (x /))l/2 ( 4 °) 

It should be noted that the profile 7(1) depends on /? and b, but since j(x) S> 1 except 
in the vicinity x ~ 1, one can easily show via integration for practically any p and b values 
that to a very good approximation the integral equals (xt — 1) and thus: 



rx'=x T j / 

T osc = (2.55 • W~ 27 s) p- 1 / = (2.55 • l(T 27 s) ^(xt - 1) (41) 

Jx'=i (1 - 7 )) v 

where is the maximum (terminal) value of x during an oscillation (Fig. 5). 

Figure 5 shows the momentum profiles obtained for p = 1 (thus applicable to both cases 
A and B) and various values of b (= £(m p /2m ) 4 /16(3). Thus increasing b corresponds 
to decreasing m and the value b = 8/15 = 0.5333 corresponds, using £ exp (Table 1), to 
m = 0.3675 eV/c 2 . The p(x) profile exhibits a maximum for all b > 6 m j n = 0.136244 but 
particle confinement occurs only for b > b c = 0.32 = 8/25, i.e. for m < 0.418 eV/c 2 . 
This is shown more clearly in Figure 6a, based on Fig. 5, which shows the dependence of 
x m (distance of maximum momentum), of xt (terminal, i.e. maximum distance during an 
oscillation) and of x p (distance of p(x p ) = 1) on b and thus m . 



m / eV/c 2 m / eV/c 2 

0.5 0.4 0.3 0.5 0.4 0.3 




Figure 6: Effect of b or of rest mass m Q on x p (p(x p ) = 1), x m (—FG(x m ) = Fc(x m )) and 
xt (terminal distance during oscillations) for p = 1 (a). Effect of b or of rest mass m on the 
period of the oscillator. Comparison of the oscillator period, when xt corresponds to the 
proton Compton length A c , with the lifetime of A baryons [22| (b). 



Figure 6b, based on (f4T)l) . shows the dependence of the oscillation period T osc on b, thus 




Figure 7: Effect of x p (distance where p(x p ) = 1, i.e. R = m p c 2 ) and corresponding 
Coulombic potential energy V c (x p ) for p = 1 on b(= (eGm 2 /4e 2 )(m p /2m ) 4 ) and on 7 P (= 
m p /2m ) showing, by the two circled points, the b and 7 P values corresponding to m = 
(l/4)(a/27r) 3 m p = 0.3765 eV/c 2 . The choice of any m and thus b(x p ) leads via this figure 
and £ = 46/ 7 4 to f = (2/15)(a/2vr) 12 (eq. (pjl ). 



m . Interestingly T osc in the vicinity of the b c = 0.32 = 8/25 limit (i.e. at the b value of 
0.320816 corresponding to xr = r X c ) equals 8.8 • 10~ 24 s, which is very close to the lifetime 



(5.6 • 10~ 24 s)of A baryons [22]] 



5.1 Maximum and critical particle rest mass 

The x p vs b curve of Figure 6a for p = 1 forms the basis, upon axis rotation, of figure 7 by 
establishing the curve b(x p ). Before presenting this figure (Fig. 7) in detail, it is useful to 
make some observations. 

In view of the definition of b (Eq. 20), the existence of 6 m ; n , as a limiting b value for 
momentum maximization ((— FG(x m ) = Fc(x m )) ()32j) . implies directly the existence of a 
minimum j p value necessary for momentum maximization, i.e. from (|20p for Q = 1/4: 

7 P ,mm = (46 mi „/£) 1/4 ; t = £ -^r L (42) 

Since from (|7|) it is m Q = m p /2j p , it therefore follows, as already known from Fig. 6a, 
that there exists a maximum particle, mass, m 0imax , above which the momentum exhibits 
no local maximum and thus particle confinement is not possible. 

The value of m max is readily obtained from: 




Figure 8: Schematic for p = 1 of the location and classical radii A 9 (= 3.06 • 10~ 3 fm) of the 
two particles at x p when b = 0.14028 (thus from Fig. 7, x p = 16) and when b = 8/15 (thus 
x p = 6.4378). In the latter case it is j p = 2(a/27r)~ 3 . 



m r , 



m n 



(43) 



" 2 7p , min 2(46 min /e) 1 /4 

Similarly the existence of b c = 8/25 as a limiting value for particle confinement, implies 
the existence of a critical m value, denoted m 0jC , above which particle confinement is not 
possible (Fig. 6a). The value of m 0jC it readily obtained from: 



m. 



2(46 c /0 1 /4 



(44) 



6 Computation of G 



So far we have not used the experimental G value, which we denote G, 



cxp 



6.674 • 



10~ n m 3 /^- 1 ^ 2 [13; |20j; [2l|. Using this value it is £ exp 
1)- 

Substituting in (|4"4l) one finds: 



0.1343595(a/2vr) 12 (Table 



2m . 



in,, 



Ccxp 

467 



1/4 



0.1343595(a/2vr) 12 
4(0.32) 



1/4 



0.5692(a/2vr) ; 



(45) 



m 0iC = 0.2846(a/2vr) 3 m p = 0.4183 eV/c 2 



(46) 



7p, c = 1.7568(a/2vr)^ ; b = b c = 0.32 = 8/25 (47) 

The computed m 0)C value lies very close to the estimated heaviest neutrino mass of 
0.4 eV/c [22| and this is quite encouraging, although neutrinos are neutral. On the other 



hand (Fig. 7) the x p value, denoted x PtC , corresponding to b c is 8.787 and the Coulombic 
energy V c is 106.78 MeV, quite low in comparison to the QCD transition energy of Eqcd = 
151(6) MeV. One then notes that since from (|20p it is: 

7 P (x p ) _ fb(x p )\ 1/A 



7p(^p,c) V b c 

for jp(xp) = 2(a/27r)~ 3 it is b(x p ) = 0.5333 = 8/15 which corresponds to x p = 6.438, 
thus (top axis in Fig. 7) Vc = 145.74 MeV, which lies within the current uncertainty limits 



of Eqcd 17]. It is also quite close to the rest energy (139.6 MeV) of ir + mesons [22j; |23j; |24j 
which are stable particles and follow electromagnetic decay 

These observations suggest a clear strategy for selecting the appropriate m value below 

The energy of Eqcd(= 151(6) MeV) of gluons at the transition temperature of QCD 
; [3] defines the minimum allowed distance, x p , of the two particles in view of their 
Coulombic repulsion, i.e. 

Vc,o/x p = Eqcd ; x p = pm p c 2 / E QC d (49) 

which for p = 1 and Eqcd = 145.74 MeV gives x p = 6.438. This x p value corresponds 
to b = 0.53333 = 8/15 (Fig. 7). In view of x = 8p(r/X q ) (eq. (5)) , this for p = 1 gives 
(r/Xq) = 0.8047 (Fig. 8), which implies, as is reasonable to expect since a bound state 
is formed, a significant overlap between the classical radii of the two particles (r p /X q = 2 
corresponds to no overlap, Fig. 8). 

We thus select m Q = (l/4)(a/2n) s m p = 0.3675 eV/c 2 , which is also very close to the 



heaviest neutrino mass of 0.4 eV/c 2 [22j|. Since 2m ~/ p = m p , this implies 7 P = 2(a/2n) 3 as 



can also be seen in Fig. 7, where the ^ p {x p ) curve has been constructed from this equality 
and ([351) . 

Thus from the definition of b(= £7„/4) one obtains: 



eGmtp _ 46 _ 4(8/15) 2 



(a/2Tr) u (50) 



e 2 7 4 2 4 (a/27r)- 12 15 

It is worth noting that in view of (|48p and £ = Ab/^f p , the same £ value is obtained by 
choosing any point x p in Figure 7. Therefore from ([501) : 



G 



2 p 2 
-(a/27r) 12 6 
15 w ; 



emi 



(2/15)(6.02370395)10- 36 (0.82464)10 : 



-,26 



6.62318 • lQ" n m 6 /kg 



-V 2 



(51) 



mental value of 6.674 • 10 11 m 3 kg 1 




which is the expression derived in |ll|; Il2| and is in excellent agreement with the experi- 



20; 121 • 



The agreement becomes quantitative if one accounts for the fact that practically all 
measurement of G are carried out using torsion balances which utilize metal rods con 
sisting of heavy elements such as Fe or Mo [19 



20 



the average nucleon mass in 56 Fe, i.e. m 



211 ] rather than protons. Thus using 
1.666157 • 10" 27 kg= 934.547 MeV/c 2 or 



that of 9S Mo(934.647 MeV/c z ) rather than the proton mass m p = 1.67262171 • 10~ 27 kg 
= 938.272 MeV/c 2 , one obtains: 



G cxp = G(m p /m Pim ) 2 = 6.6760 • 10~ n m 3 kg- 1 S - 2 (Fe) 

= 6.6746 • 10 -11 rn^kg^s' 2 (Mo) (52) 

which is in quantitative agreement with experiment (Fig. 9a). 
Conversely when solving (|51|) for m Ptm one obtains: 

2 \ 1/2 / 2 x V2 

(a/2vr) 6 — = 1.666 • 10" 27 kg = 934.6 MeV (53) 



15/ v ' ' \eG) 

which is also in excellent agreement with experiment (Fig. 9b) particularly for those 
elements (Fe, Mo, Cu) commonly used in measuring G exp with metal rod torsion balances 
0; 0; 21]. As Fig. 9b shows ([52]) predicts a small (~ 0.01%) variation in the G value 



measured with different metal rod torsion balance elements and such a small variation has 



indeed been observed for years [19l;l2ll]. A prediction of (|52|) . (|53|) and Figure 9b is that if one 
could somehow measure G~exp using protons or U2 instead of metal rods, then the measured 
value would be 6.62318 ■ 10~ n m 3 kg~ 1 s~ 2 , i.e. 0.8% smaller than the experimental value, 
but this is obviously a very difficult experiment. 

6.1 Maximum m value 

From (|4"3|) and the computed value of ^ (eq. 50) one can compute the value of m 0jmax . Thus: 

m , max = -3?- = m P > = 0.351(Q/2vr) 3 m p = 0.516 eV/c 2 (54) 

Similarly one can use the critical value of b for particle confinement, i.e. b c = 0.32 
to compute the corresponding critical particle mass m 0)C above which particle confinement 
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Figure 9: Predicted and measured gravitational constant and average nucleon mass, (a) 
Comparison of (|5ip with the time evolution of the experimental gravitational constant value, 
and with some of the most recent experimental values [19,21]; Pre-1997 values data from 
Table 2 in Ref. 21. Post-1997 values from Table X in Ref. 19. Shaded areas shows the 
CODATA recommended [19] value of (6.6742 ± 0.001) x 10~ n m 3 kg- 1 s- 2 (b) Comparison 
of the average nucleon mass in various nuclei with the one predicted by equation (1521) . 
The right ordinate shows the value of G ex p expected to be measured with different torsion 
balance rod materials. 



does not occur: 



m. 



m r . 



0.2814(a/27r) 3 m p = 0.4183 eV/c" 



(55) 



2 7p , c 2(4b c /0 lfi 

As already noted, the value of m a corresponding to 7 P = 2(a/27r)~ 3 is 0.3675 eV/c 2 , i.e. 
is smaller than both 

^o,max cind TYi c . All thes e va lues are in excellent agreement with the 
estimated heaviest neutrino mass of 0.4eV/c |22|. 



In brief, although all m values below m 0jC lead to the formation of confined states 
(Figures 5 and 6) and via (|50|) and ([52]) and Figure 7 lead to the (wp value of 6.674 • 
10~ n rr?kg~^ s~ 2 in quantitative agreement with experiment 19; l20i; 121], the selection of 
m = (l/4)(a/2-7r) 3 m p also leads to: 

1. A TUo value (0.3675 eV/c 2 ) quite close to the estimated heaviest neutrino mass of 
0.4eV/c 2 0. 

2. A V c energy (145.2 MeV) in good agreement with the QCD transition energy of 
151(6) MeV and also with the energy 139.6 MeV, of vr+-mesons Q. 



7 Consistency with quantum mechanics 

The above analysis was based entirely on special relativity and on the laws of Coulomb and 
Newton. Thus for completeness the results should be compared with quantum mechanics. 

As a first step it is interesting to note that the confinement distance xtt of the two 
particles lies below or close to the proton Compton length, A c , and thus below or close to 
their de Broglie wavelength, which is something one normally would expect to find from 
quantum mechanics. Thus for the particle with m Q = (l/4)(a/27r) 3 m it is b = 8/15 and 
7 P = 2(a/27r) -3 (Fig. 7) and thus the particle momentum p(x p ) at x = x v is given by: 

p(x p ) = J p m c = 2(a/27r)~ 3 (l/4)(a/27r) 3 m p c = (l/2)m p c (56) 
Thus the particle de Broglie wavelength, A p , is given by 

A »=(I72^ = 2A " (57) 

i.e., it coincides with the Compton wavelength of a particle with mass m p /2, in consis- 
tency with the uncertainty principle. 

The agreement between the maximum separation xxr and A c becomes quantitative for 
b = 0.320816 w 8/25 (Fig. 10 which is based on Fig. 6b and eq. (|4*T|) ). Very interestingly 
in this case it is T osc = 8.81294 • 10 -24 s thus / = 1.1347 • 10 23 s" 1 and u = 2vr/ = 
7.1295 • lO 23 ^ 1 , therefore hf = huj = 7.5186 • 10~ n J = 469.274 MeV = (l/2)m p c 2 (Fig. 
6b and Fig. 10), which is the total energy, E, of each particle in the confined state. Thus 
in this case two equations follow: 




Figure 10: Effect of xt and oscillation amplitude on the period, frequency and energy 
of the oscillator. 



h = (c/f){m p c/2) = ArAp (58) 

where Ar and Ap express the uncertainty in distance and momentum due to the oscil- 
lations and: 

E = (l/2)m p c 2 = hf (59) 



The first equation (|58p shows the exact conformity of the solution, /, obtained from 
special relativity with Heisenberg's uncertainty principle and the de Broglie equation, while 
the second equation (|59p shows that, interestingly, the total energy of the oscillating particle 
is expressed in terms of / by the Planck equation for the energy of a photon. 

It is also interesting to note that the energy, E, corresponding to b = 8/15, thus T osc = 
2.876- 1(T 26 s (Fig. 10) is 143.8 GeV which lies within the Standard Model predictions about 
the mass (115 — 180 GeV/c 2 ) of the Higgs boson 



25 



26J, but this may be a coincidence. 
Finally it is worth noting that in view of the very small value of m c 2 (~ 0.4 eV) the 



m c 2 term in (fSj) i.e. 



E(x, t) = yjm\& + p 2 (x, t) ■ c 2 (60) 

is negligible and thus it is E(x, t) pa p(x, t) ■ c which implies that the time dependent 
Klein-Gordon equation reduces to the time-dependent Schrodinger equation: 



*£*< r -'> = (-lb5? +vw '•<'•'> (61) 



which in turn reduces, due to the very small m value, to: 



+ 9 t, \ & d 2 ^(r,t) , . 



which is the free-particle wave equation with solutions of the form 

V(x,t) ~ e i ( kr - ut ') (63) 

where k = uj/c ~ 2.37 • 10 15 m -1 , i.e. k = 2.37 fm _1 when using the above determined 
<jj value. Thus in summary one may conclude that the creation of the bound oscillatory 
state by the two particles, obtained from special relativity, is also consistent with quantum 
mechanics. 



8 Gravitational confinement in circular trajectories 

The two light particles of mass m Q forming bound hadron states need not be charged if their 
initial kinetic energy is sufficiently high. This point is worth analyzing, since neutrinos, 
e.g. electron neutrinos, have a charge radius < r 2 >< 3.32 • 10 _36 m 2 and < \r Ve \ >< 
1.83 • 10 _3 /m [27,28]) but no net charge. 

Thus if the two particles of rest mass m have velocities v and -v relative to the laboratory 
observer and are moving initially in opposite directions on two parallel paths sufficiently 
close to each other, then , to a good approximation, the gravitational force judged by the 
laboratory observer is again: 

Gm 2 i & 

F G = ^f- (64) 

where r is again the distance and this force acting between the two particles becomes 
the centripetal force for their cyclic motion, i.e. [10]: 

Gm 2 -f 6 _ -fm c 2 _ GffioT 5 _ Gm n rf 

r 2 ~ r/2 ; r/Z ~ c 2 ~ 2c 2 ^ 

where m n is the neutron mass and in the last equality we have set 7 = 7 P (which 
guaranties that \x = 1 and thus the confined mass is m n ) and have used m n = 2j p m 
similarly to ([7]). The problem described by (|65p is very similar to the Bohr treatment of 
the H atom with the gravitational attraction replacing the Coulombic attraction. 

In the linear motion problem we have chosen the parton mass from m = (l/4n^ t )(a/2-7r) 3 m p 
with n m = 1. In the circular motion we choose n m = 2, thus m a = (l/16)(a/27r) 3 m n and 
thus 7 P = 8(a/27r) -3 . 



As in the Bohr treatment of the H atom we assume that the angular momentum is 
quantized, i.e. 



7 p m c(r/2) = nh (66) 
and using to express r/2 and accounting for 7 p = 8(a/27r) -3 one obtains: 

G= McW2< thus for n = l G = 6.6293 • 1CT 11 rr^kg^s^ 2 (67) 

which very interestingly, differs less than 0.1% from the value computed in (|5ip . which 
involves m p instead of m n , and less than 0.4% from the experimental G value of 6.6742 • 
10 -11 m^kg^s- 2 [19,20,21]. 

The mass m of each one of the rotating particles is m = (l/16)(a/27r) 3 m p = 0.0919 
eV/c 2 which is within the range of neutrino masses [22]. The distance r of the two particles, 
i.e. the diameter of rotation, is obtained from (|65|) . i-e. 

Gm„8 4 (a/2vr)- 12 
r = — V '- = 0.844 fm 68 

c z 

which is in excellent agreement with the estimated proton or neutron diameter A c = 
1.32 fm. The period of the rotation is irr/c = 8.838 • 10 -24 s, i.e. almost identical with the 
period 8.813 • 10~ 24 s computed for the linear motion with r iy = A c (Fig. 6b) and very 
close to the lifetime of A baryons (5.6- 10 -24 s, Fig. 6b). Thus one may safely conclude that 
gravity can confine both charged and neutral particles of the mass range of neutrinos to 
form hadrons. The same analysis presented in this section also applies to pairs of charged 
and uncharged light particles (e.g. pairs with charges e and zero) and this eliminates the 
necessity of assuming fractional charges Qx, Q2 (e.g. 1/2 or ±2/3 and ±1/3) as in the linear 
motion case to form charged hadrons. 



9 Discussion and conclusions 

The analysis of the one-dimensional motion of two charged particles under the influence of 
their Coulombic repulsion and gravitational attraction using the relativistic equation of mo- 
tion has shown that when the particles are sufficiently light (m G < 0.4183 eV/c 2 , in the mass 
range of neutrinos) then gravity is sufficient to create a stable electrostatic-gravitational os- 
cillator, i.e. a bound state, having the mass of a proton. This is also consistent with the 
de Broglie wavelength expression and with the Klein-Gordon and Schrodinger equations of 
quantum mechanics. 

The transformation of Coulombic energy into kinetic energy of the two particles and thus 
into rest energy R (and rest mass m = R/c 2 ) of the confined state, provides a straightforward 
mechanism of "hadronization" [22], i.e. of generation of rest mass when a bound state is 
formed by two fast particles. 



Also when the initial kinetic energy of two neutral particles (neutrinos with mass 0.0919 
eV/c 2 ) is sufficiently high (m p c 2 ), gravity alone suffices to confine these particles in stable 
circular (or elliptical) trajectories corresponding to neutral hadrons, e.g. neutrons. The 
relevance of the present simple analysis to the actual problem of hadron formation via 
the condensation of quark-gluon plasma [13, 16-18] appears worth further investigation for 
several reasons: 



1. It is interesting and reassuring that the computed maximum particle mass m 0)C ~ 
0.4183 eV/c 2 necessary for confinement lies in the range of the heaviest neutrinos 
(~ 0.4 eV/c 2 ). This shows that such small particles as the ones used in the present 



221. It is 



analysis, actually exist and in fact have velocities near the speed of light 
also worth reminding that during the decay of muons two, rather than one, neutrinos 
are produced [22]. 



2. The two fast moving (or rotating) particles appear to have the basic properties of 



gluons, i.e. they induce the strong force and they are practically massless [22j; |23j; 



24 



29i ; [30J. Thus the formation of neutral hadrons such as neutrons, appears to 
follow logically from the present analysis. For the case of charged hadrons it follows 
that if one may view gluons as charged neutrinos, then the creation and stability of 
protons could also be readily rationalized on the basis of the present analysis. As 
already noted, in the case of circular trajectories the necessity of assuming integer or 
fractional charges is eliminated. 

The massive particles formed during each oscillation (as (jl{x) oscillates and thus 
the rest mass of the system oscillates between 2m Q (~ 0.8eV/c 2 ) and 2ym (> 
932MeV/c 2 )) appear to have some of the basic properties of sea quarks [22J, i.e. 
they are produced and annihilated as virtual particles during each oscillation, exactly 
as envisioned in the standard model [22| . The maximum period of these oscillations 



(~ 8.7 • 10 24 s) lies very close to the lifetime (~ 5.6 • 10 24 s) of A baryons [22 

In this direction it is worth pointing out that the present analysis of the motion of 
particle B is directly applicable, with very minor modifications, to the motion of the 
same particle in presence of two other particles in a symmetric arrangement (Fig. 11), 
and this geometry, to which the present results are directly applicable, involving three 
particles, is very close to the standard model picture of three quarks produced and 
annihilated continuously inside a proton 221 ] . 



4. The Coulombic and the negative of the gravitational potential energy at the con- 
finement point x p are equal to 145.7 MeV, very close to the energy 151(6) MeV of 
marks and gluons at the condensation point of quark-gluon plasma to form hadrons 



qu£ 





; 0; [l?]; Q] and also close to the energy, 139.6 MeV, of 7r + -mesons 




Figure 11: Geometry of the one dimensional motion of particle B involving two other 
Coulombically repelling charged particles Al and A2 replacing particle A in Figure 1 in a 
symmetric arrangement. The two reference frames S and S' remain the same. 



5. The limiting (for r — > oo) escape rest mass R/c 2 for b < 8/25 (e.g. b = 2/9, Fig. 3a) 
is very near the rest masses of all the (1/2) spin hadrons, i.e. the £, £ + , S° and X - 
baryons [22j] (Fig. 3a). 



6. The rest mass R/c 2 for b > 8/25 (e.g. b = 8/15, Fig. 3b) during the oscillations 
between x p and xt covers the range of the rest masses of all the (3/2) spin baryons, 
i.e. the A, £*, S* and 0™ baryons (Fig. 3b). 

7. The results are consistent with quantum mechanics, and when v xt = A c , then the 
energy Twj equals (l/2)m p c 2 (Figs. 6b and 10). 

8. Finally the analysis provides a straightforward explanation about why individual 
quarks or gluons cannot be separated from hadrons or mesons and studied indepen- 
dently. Although the analysis shows that, at first surprisingly, the strong force can be 
viewed as the relativistic gravitational force, the emerging picture is very similar to 
that of the standard model regarding the strong force and the composition of hadrons 



22 



23 



29i ; 1301 ] as well as their behaviour in elastic and inelastic scattering experi- 

of the bag 
30]. 



ments which has led to the concepts of partons connected by strings 
model 2J] and of the massive quarks and practically massless gluons 



22 



23 



29 



All these concepts are consistent with the present analysis at different phases of the 
oscillation. 



Regardless of the exact direct or indirect relevance of the present simple analysis to the 
above important physical phenomena, particles, and concepts it is certain that gravitational 



forces suffice to create stable electrostatic-gravitational oscillator states and this leads in a 
straightforward manner to simple formulae (|51l to[53landl67 |) for the gravitational constant 
and for the proton mass in terms of the other physical constants which are in quantitative 
agreement with experiment. 

It is also reasonable to expect that the gravitational forces exerted between the fast 
moving particle constituents (partons or quarks) of neighboring hadrons at fm distances 
can also be quite important and thus can lead to binding energies per nucleon of the order 
of arripC ~ 7 MeV [22| and the formation of nuclei. Preliminary work has shown 

that indeed the binding energies of some light nuclei ( 2 H and A He) can thus be computed 
using (|5ip with good accuracy and this point also appears to deserve further investigation. 
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Table 1. CODATA recommended values of e, e, m p and G and 
comparison of the value of £(= eGm 2 /e 2 ), (= £ exp ) computed 
from them and from equations (|50|) . (|51|) and (1 521) 

e = 1.6021765 x 10" 19 C 

s = AttEo = 1.112649 x 10~ 10 C 2 /Nm 2 

h = 6.6260693 x 10" 34 Js 

c = 2.997925 x 10 8 ms" 1 

m n = 939.565 MeV/c 2 = 1.67492728 x 10~ 27 kg 

m p = 938.2723 MeV/c 2 = 1.67262171 x 10~ 27 kg 

Gexp = 6.6742 x 10" 11 r^kg^s' 2 

^cx P = 8.093421 x 10~ 37 = 0.1343595 x (a/2vr) 12 

£ = 8.031608 x 10~ 37 = 0.13333 x (a/2vr) 12 

G = 6.62318 x 10- 11 m 3 ^"^- 2 = 0.13333(a/2vr) 12 (e 2 /em 2 ) 

G = 6.674 x lO^m 3 ^ -1 * -2 = 0.13333(a/27r) 12 (e 2 /em 2 m ) 

m P)m in 56 Fe = 934.547 MeV/c 2 and in 98 Mo = 934.647 MeV/c 2 

m PtTn from eq. ([53j) : m PtTn = 934.6 MeV/c 2 

a = 1/137.035 

(a/271-) = e 2 /ec/i = 1.16141822 x 10~ 3 
(a/2?r) 3 = 1.5666281 x 10" 9 
(a/2?r) 12 = 6.02370395 x 10~ 36 
(a/2ir)- 12 = 1.6601081 x 10 35 



Table 2. List of symbols 

b: ( e Gm2/16Qe 2 ) 7p 4 = e(7 P 4 /16Q) 
rug : particle longitudinal mass, 7 3 m 
m : particle rest mass 
m n : neutron mass 
m p : proton mass 

fnp,m : average nucleon mass in a nucleus 

n,n m : positive integers 

p : momentum 

q± , (72 : particle charges 

Qi, Q2 : <?i/ e > 92/e 

Q: Q1Q2 

x : dimensionless distance r/r Q 

x m : value of x at maximum particle momentum 

x p : value of x at /x(x) = 1 

xt '■ terminal x value during oscillations 

y(x): -V G (x)/V c ,o 

Greek symbols 

a : e 2 /ech 

(a/2n) : e 2 /ech 

(5: v/c 

7: (l-^/c 2 )- 1 ^ 

7 P : value of 7(x) at x = x p 

A c = hjvTLpC 

X q = e 2 /e(m p /2)c 

fj,(x) : dimensionless mass = j(x)/^ p = R{x)/m p 
£: (eGm 2 p /e 2 ) 

P ■ Vc,o/"i P c 2 = Qe 2 /er m p c 2 



Table 3. Summary of key mathematical equations describing the energy balance and 
the particle equation of motion 





Case A (2K = (1 - p)m p c 2 ) 


Case B (K Q = 0) 


Energy balance and 
corresponding 
gravitational potential 
energy profile 

tfW V c ,o 


dy(x) _ ^(l-p/x+ P y(x)) 6 i^oj) 


= 46p 6(l-V^+3/W) 6 (J23J) 






BC: 2/(1) = 


H{x) = l-p{{l/x)-y{x))m 


//(a:) =/3(l-l/x + y(x)) ([25]) 


Equation of Motion 


| = F c (x) [1 - 4^ 6 (x)] = (2p 2 m p c 2 /QX q x 2 ) [l - 46/i 6 (x)] flE}, 42D 


T osc « (A g /4pc)(x T - 1) (0DD, (0U) 



